From E 2 n to other modes : Effects of pressure on electron-phonon interaction in MgB 2 
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We study the effects of pressure on the electron-phonon interaction in MgB^ using density- 
functional-based methods. Our results show that the superconductivity in MgBi vanishes by 100 
GPa, and then reappears at higher pressures. In particular, we find a superconducting transition 
temperature T c fa 2 K for fj,* = 0.1 at a pressure of 137 GPa. 



The superconductivity in MgB 2 1] is intimately con- 
nected with the holes in the B p x ( y ) and p z bands 0.0 0| 
and their coupling to the in-plane E 2g phonon mode 
a h ll al . Several experiments on MgB 2 under pressure 
lfj,|ri| show a reduction in superconducting transition 
temperature T c at a rate of -1.1 K GPa , which can be 
understood in terms of a continuous weakening of E 2g 
phonon mode coupling. Based on this picture, it is es- 
timated that the superconductivity in MgB 2 will vanish 
by 90-100 GPa [B]. 

Recent experiments on elemental Li 0,0] and Y 0] 
have shown a significant increase in their T c with increase 
in pressure. For example, the T c of Y changes from 6 mK 
Efj to 20 K as the pressure is increased from ambient 
to 115 GPa. In light of these experiments, it would be 
interesting to see what happens to the superconducting 
properties, especially T c , of MgB 2 beyond 90-100 GPa. 

In this Letter, we report on our density-functional- 
based study of the evolution of electron-phonon inter- 
action in MgB 2 as a function of pressure up to 150 GPa. 
Additionally, we have solved the isotropic Eliashberg gap 
equation to obtain T c of MgB 2 as a function of pressure. 

Based on our calculations, we find that the increase 
in pressure hardens the MgB 2 lattice and dramatically 
reduces the contribution of the E 2g phonon mode to the 
electron-phonon coupling. As a result, the superconduc- 
tivity in MgB 2 vanishes by 100 GPa, only to reappear at 
higher pressures. In particular, we find a superconduct- 
ing transition temperature T c w 2 K for //* = 0.1 at a 
pressure of 137 GPa. We also find that the B\ g phonon 
mode completely softens around A symmetry point by 
143 GPa. Before we describe our results, we briefly out- 
line the computational details of our calculations. 

We have calculated the electronic structure and the 
electron-phonon interaction of MgB 2 in P6/mmm crys- 
tal structure as a function of pressure up to 150 GPa. The 
lattice constants a and c for different pressures were taken 
from Ref. [12], which agree well in the range of pressures 
where experiments have been reported. The electronic 
structure was calculated using the full-potential linear 
muffin-tin orbital (LMTO) method [H [3] as well as 
the plane-wave pseudopotential method using PWSCF 
package The phonon dispersion, Wq„, the phonon 

linewidth, 7 q „, the phonon density of states, F(ui), and 
the Eliashberg function, a 2 F(uj), were calculated using 
the density-functional perturbation theory jijj] as imple- 



mented in the PWSCF pac kage. Finally, the isotropic 
Eliashberg gap equation |2lL Ea has been solved for 
a range of n* to obtain the corresponding T c . 

The charge self-consistent, full-potential, LMTO calcu- 
lations for MgB 2 were carried out using the local-density 
approximation for exchange-correlation of Perdew et al. 
[2J|, 3/t-energy panels, and 793 k-points in the irreducible 
wedge of the Brillouin zone (BZ). The 2p state of Mg 
was treated as a semi-core state. The basis set used 
consisted of s, p, and d orbitals at the Mg and the B 
sites, and the potential and the wave function were ex- 
panded up to Imax = 6. The muffin-tin spheres for Mg 
and B were taken to be slightly smaller than the touch- 
ing spheres in each case. These calculations are used to 
study the symmetry-resolved densities of states as well as 
the band-structure along the high symmetry directions in 
the corresponding BZ. 

For calculating the Fermi surface, the phonon disper- 
sion, the phonon density of states, and the Eliashberg 
function, we have used the plane-wave pseudopotential 
method with Vanderbilt's ultrasoft pseudopotentials with 
nonlinear core correction. The kinetic energy cutoff for 
the wave function was 35 Ry, and for the charge den- 
sity and the potential it was 200 Ry. The exchange- 
correlation potential was parametrized as suggested by 
Perdew et al. 0|. The Fermi surface was constructed 
using eigenvalues calculated on a 36 x 36 x 36 grid in the 
reciprocal space. The linear-response calculations were 
carried out on a 6 x 6 x 6 grid resulting in 28 irreducible 
q-points. A 12 x 12 x 12 grid of k-points is used for 
the BZ integrations during linear response calculations. 
The electron-phonon matrix elements are sampled more 
accurately on a 36 x 36 x 36 grid using the double-grid 
technique. To check the convergence of our linear re- 
sponse calculations, we have used an 18 x 18 x 18 and 
24 x 24 x 24 grids for carrying out the BZ integrations 
during self-consistency for points along T~A direction as 
well as for pressures above 125 GPa. The convergence of 
electron-phonon matrix elements were checked with grids 
up to 40 x 40 x 40 in the reciprocal space. 

The B p x (y) and p z densities of states (DOS) play a 
crucial role in determining the superconducting proper- 
ties of MgB 2 . Our calculated B p x ( y ) and p z DOS as a 
function of pressure are shown in Fig. where we also 
show the DOS in a 10 mRy energy window around Fermi 
energy, Ejr, in a separate panel. With increasing pres- 
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Figure 1: The B p x ( v ) (solid line) and the p z (dashed line) 
densities of states in MgBi as a function of pressure for 
(1st panel from the top), 50 (2nd panel), 100 (3rd panel), 
and 137 (4th panel ) GPa. The panels on the right show the 
corresponding densities of states in a 10 mRy energy window 
around the Fermi energy, which is shown as the vertical, long- 
dashed line in the figure. 
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Figure 2: The band-structure of MgBi along high symmetry 
directions in the Brillouin zone as a function of pressure for 
(1st panel from the top), 50 (2nd panel), 100 (3rd panel), and 
137 (4th panel) GPa. The horizontal line, passing through 
the energy zero, indicates the Fermi energy. Note that the 
numerical values of the wave vector k corresponding to the 
various symmetry points, as shown in the figure, are different 
at different pressures. 




sure both p x t y ) and p z DOS flatten and the states are 
pushed further down. Since with increasing pressure the 
lattice constant c decreases much more rapidly than a in 
MgB 2 , as a result the p x ( v ) DOS around Ep is signifi- 
cantly reduced. Between 100 and 137 GPa the p x (y) DOS 
is inside Ep with essentially no holes. However, the p z 
DOS around Ep is in tact up to 137 GPa. Such a sit- 
uation is very different from the usual band filling that 
takes place when we put Al and/or C as impurities in 
MgB 2 - These differences can be further highlighted by 
considering the changes in the band-structure along high 
symmetry directions in the BZ of MgB 2 as a function of 
pressure. 

In Fig. [21 we show the calculated band-structure of 
MgB 2 along high symmetry directions in the Brillouin 
zone. The states along T-A have B p x / y ) character while 
states at M are of B p z type. With increasing pres- 
sure, as can be seen from Fig. [2 the states along T-A 
move towards Ep while the states at M move away from 
Ep. By 137 GPa, all the B p x ( y )-derived states along 
T-A are close to being inside Ep, while the states at 
M have moved up by almost 0.17 Ry. The movement 
of states along T-A and M with increase in pressure, in 
conjunction with doping of MgB 2 , can be more effective 
in tuning its superconducting properties. 

A complete picture of the pressure-induced changes in 
the electronic structure of MgB 2 emerges through the 
changes in the Fermi surface, as shown in Fig. We find 
that the characteristic two-dimensional cylindrical sheets 
along T-A become narrower and constricted around T 
point with increase in pressure. With further increase in 
pressure, say by 100 GPa, the states manage to form only 
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Figure 3: The Fermi surface of MgB2 as a function of pressure 
in the corresponding Brillouin zone for (a) 0, (b) 50, (c) 100 
and (d) 137 GPa, calculated as described in the text. 
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Figure 4: The phonon dispersion uo of MgBi along high sym- 
metry directions in the Brillouin zone (left panel), and the 
w q „ (top-right panel) and the 7 q „ (bottom-right panel) for 
Bi g phonon mode along T-A-L, as a function of pressure for 
125 (dashed line), 137 (solid line) and 143 (dot-dashed line) 
GPa. The Uq V and the 7 q „ are in meV, the frequencies below 
zero (in the top-right panel) are imaginary, and for 143 GPa 
the 7q„ is not shown. Note that the numerical values of q 
corresponding to the various symmetry points are different at 
different pressures. 



a small portion of the tube around A point of the T-A 
direction. The circular opening of the p z -derived Fermi 
surface sheets containing {K,M) and (H,L) symmetry 
points narrows as the band moves up. At a pressure of 
137 GPa, the Fermi surface consists of essentially these 
sheets with the cylindrical tube reduced to a dot around 
A point. Thus, as long as the states on the cylindrical 
sheets play a pivotal role in sustaining superconductiv- 
ity in MgB 2 , their absence beyond 100 GPa, as shown 
in Fig. 03 ensures the vanishing of superconductivity in 
MgB 2 around 100 GPa and above. Unless, of course, 
with pressures above 100 GPa the electron-phonon cou- 
pling provided by the remaining states is strong enough 
to restart the superconductivity in MgB 2 . 

To be able to go beyond the qualitative assessments 
that we have made so far about the pressure-induced 
changes in the electron-phonon interaction in MgB 2 us- 
ing the changes in the electronic structure as the back- 
drop, we have calculated the phonon dispersion, the 
phonon linewidth, the phonon density of states and the 
electron-phonon interaction as a function of pressure up 
to 150 GPa. We find that in MgB 2 at zero pressure the 
frequency associated with E2 g phonon mode at T is 67 
meV, which goes up to 148 meV at 137 GPa. Similarly, 
the frequency of the E2 U phonon mode at A point in- 
creases from 60 meV at zero pressure to 139 meV at 137 
GPa. In contrast, the frequency of the B\ g phonon mode 
at A point, involving out-of-plane motion of B atoms, 
reduces from 63 meV at 100 GPa to 34 meV at 125 GPa. 
As shown in Fig. Q] the softening of B\ g mode around A 
point continues with further increase in pressure and by 
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Figure 5: The phonon density of states F(lo) (top 2 pan- 
els) and the Eliashberg function a 2 F(u>) (bottom 2 panels) of 
MgB2 as a function of pressure. The F(uj) and a 2 F(u>) for 
(solid line) and 50 (dashed line) GPa are shown in the 1st 
panel and the 3rd panel, respectively. The F(u>) and a 2 F(u>) 
for 100 (dotted line), 125 (dashed line) and 137 (solid line) 
GPa are shown in the 2nd panel and the 4th panel, respec- 
tively. 



143 GPa it has become unstable. The softening of B lg 
mode around A point is accompanied by an enhanced 
7qi/, as can be seen from Fig. 0J Note, however, that 
the dynamical instability of B\ g mode in the present cal- 
culation does not necessarily imply a physical instability 
of the lattice. Under these conditions, the anharmonic 
effects become important. 

Our calculated F(uj) and a 2 F(uj) of MgB 2 for pres- 
sures up to 137 GPa are shown in Fig. For MgB 2 
at zero pressure, our results for F{to) compare well with 
previous calculations HJ 0- The decrease in the 
strength of the electron-phonon coupling in MgB 2 with 
increase in pressure can be clearly seen in Fig. Q3 The 
average electron-phonon coupling constant A, shown in 
Fig. El decreases from 0.68 to 0.28 as the pressure is in- 
creased from zero to 100 GPa, which is consistent with 
the expected loss of superconductivity in MgB 2 by this 
pressure. By 100 GPa the partial X qv of the E 2g phonon 
mode at T and A points are reduced by a factor of 3 and 
20, respectively. However, the increase in A to 0.38 by 
137 GPa with the possibility of higher values with fur- 
ther increase in pressure, as shown in Fig. suggests 
the reemergence of superconductivity in MgB 2 at higher 
pressures without the dominance of E 2g mode. 

Finally, the change in T c with pressure has been ob- 
tained by solving the isotropic gap equation (12, |3] using 
the calculated Eliashberg function a 2 F(ui) at each pres- 
sure. The results of these calculations are shown in Fig. 
HO for a range of /i*. Our results, shown in Fig. |£j predict 
the reemergence of superconductivity in MgB 2 at pres- 
sures above 100 GPa. In particular, we find that at 137 
GPa MgB 2 will superconduct with T c w 2 K for /i* =0.1. 
Since the electron-phonon interaction in MgB 2 at zero 
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Figure 6: The calculated superconducting transition temper- 
ature T c in MgB-2 as a function of pressure P for equal to 
0.1 (solid circles connected with solid line), 0.15 (solid squares 
connected with dotted line) and 0.2 (solid triangles connected 
with dashed line), respectively. The inset shows the variation 
of A with pressure. The A and T c at 143 and 150 GPa are 
obtained by ignoring the coupling of imaginary frequencies. 



pressure is known to be very anisotropic, an accurate de- 
termination of T c requires the solution of the anisotropic 
gap equation. However, we have seen that the B p x (y)- 
dependent electron-phonon coupling decreases dramati- 
cally with increasing pressure, making the use of isotropic 
gap equation for the calculation of T c more reliable at 
higher pressures. 

In conclusion, we have studied the effects of pressure on 
the electron-phonon interaction in MgB 2 using density- 
functional-based methods. The increase in pressure hard- 
ens the lattice and dramatically reduces the contribution 
of the E 2g phonon mode to the electron-phonon coupling. 
As a result, the superconductivity in MgB 2 vanishes by 
100 GPa, only to reappear at higher pressures. In partic- 
ular, we find a superconducting transition temperature 
T c k 2 K for /j* = 0.1 at a pressure of 137 GPa. 
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